CHARGE FLUCTUATIONS IN A QUANTUM POINT CONTACT ATTACHED TO 

A SUPERCONDUCTING LEAD 
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We show how to calculate the charge noise spectrum in a normal mesoscopic conductor, which 
is capacitively coupled to a macroscopic gate, when this conductor is attached to L normal leads 
and M superconducting leads, the only restriction being that the superconducting leads must be 
at the same chemical potential. We then proceed to examine results for a quantum point contact 
(QPC) in a normal lead connecting to a superconductor. Of interest is the fluctuating current in 
a gate capacitively coupled to a QPC. The results are compared with the case when all leads are 
normal. We find a doubling of the equilibrium charge fluctuations and a large enhancement (> 2) 
in the current noise spectrum to first order in \eV\, when a channel in the QPC is opening. 
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Theory and experimental measurements of the electri- 
cal transport properties of mesoscopic conductors, which 
either contain superconducting regions or are attached to 
superconducting leads have generated great interest. Of 
particular interest are the noise properties of normal su- 
perconducting interfaces. While many aspects of the low 
frequency-current noise spectra jlj have been understood 
by generalizing the scattering approach of normal 
conductors, the fluctuations of the charge have remained 
unexplored. We generalize the work done for normal 
mesoscopic conductors to systems which contain 

superconducting leads. The only restrictions are that 
there must be one or more normal leads, the supercon- 
ducting leads must all be at the same chemical potential 
and any gates in the structure only see normal regions of 
the conductor. Having developed this technique we shall 
consider one particular example shown in Fig. |l|, a Quan- 
tum Point Contact (QPC) attached to one normal lead, 
one superconducting lead and capacitively coupled, via 
the Coulomb interaction, to a macroscopic gate. We are 
interested in the charge fluctuations in the hybrid struc- 
ture which can be measured by observing the current 
fluctuations at the gate. It turns out that at equilibrium 
the current fluctuations are determined by the RC -time 
constant. Thus our primary aim is to find the charge 
relaxation resistance R and the capacitance C for hybrid 




structures. In the presence of transport, a novel resis- 
tance Ry, which reflects the shot noise ||-| of the hybrid 
structure, appears. 

Consider a mesoscopic conductor with L normal ideal 
leads and M superconducting leads. The conductor can 
be described by a scattering matrix with elements s 7q( 5 (3 
which relate particle and hole amplitudes (a, (3 = p,h) 
incoming at contact 6 to the outgoing amplitudes at con- 
tact 7. Over a region Q the mesoscopic conductor is 
capacitively coupled to a gate. The current fluctuations 
induced into the gate can be found from the charge fluc- 
tuations of the mesoscopic conductor MM . The investi- 
gation of the charge fluctuations starts from the analysis 
of the bare charge fluctuations which are then screened 
to find the true charge fluctuations For a normal 

conductor the bare charge fluctuations are determined by 
a density of states matrix which has dimensions L * L. 
For hybrid structure the matrix has dimensions 8 * L * L: 
Any scattering channel, independently on whether the in- 
coming or outgoing channels are electron or hole like, can 
contribute to the electron density at a point r inside the 
conductor. Thus to express the density in terms of the 
scattering matrix requires a conceptual trick: It is useful 
to imagine that there are two electrostatic potentials U p 
and U h which act each separately on the electrons and 
holes [[To) . The matrix which governs the charge fluctua- 
tions can then be found by testing the scattering matrix 
elements with regard to small variations of the electron 
and hole potentials. This procedure gives the following 
density of states elements [ IO) 
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FIG. 1. Quantum point contact, attached to one normal 
and one superconducting lead, capacitively coupled to a gate. 



where the labels v : 7 and 8 denote contacts (1....L), a, (3, 
A and rj denote the electron/hole degrees of freedom (p/h) 
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and q p 



The functional derivatives are taken 



at the equilibrium electrostatic potential U p — U h — U eq . 
To give an example, Ni p2p {lh, r) is the hole density asso- 
ciated with two electron current amplitudes incident from 
contacts 1 and 2 and a reflected outgoing hole amplitude 
in contact 1. With the help of these basic expressions we 
can now find both the average density of states as well as 
the fluctuations. The hole density of states of a region 
of the conductor is N h = N h {p) + N h (h) where 



N h (a)=J2 J d 3 rTr[< 7 > A ,r)] 

1/7 A ^ 



(2) 



and the particle density of states is N p = N p (p) + N p (h) 
with 



(3) 



The trace is over open quantum channels. N a {(3) is the 
injectivity of particles (holes) (j3 = p(h)), from all con- 
tacts into the conductor, given a change in the particle 
(hole) potential (a = p(h)). 

The fluctuations of the bare charge in a region of 
interest can be found from the charge operator eJ\f given 
by 

3- V u\ Ja 

x l' W ts f 3 ("Ai r;E,E + huj)a Sfi (E + hu>) (4) 

where the zero-frequency limit of ACf s (v\, r; E, E + Tiui) 
is given by Eq. (|l|). In Eq. (|J), (E) creates an incom- 
ing electron/hole (a = p/h) in lead 7. The true charge 
fluctuations must be obtained by taking into account the 
Coulomb interaction and below we show how to obtain 
the true charge fluctuations from the fluctuations of the 
bare charges. 

Given equations |l]-|]) our next step is to consider the 
system shown in Fig. 1 and to calculate the charge fluc- 
tuations in the mesoscopic region and hence deduce the 
current fluctuations in the gate. To do this it helps to 
first consider an approach which allows us to correctly 
calculate the charge SQ in the mesoscopic region when 
a voltage is applied across the system. The charge ac- 
cumulated underneath the gate in the region Q can be 
expressed in two ways |t]-|9| : If we describe the Coulomb 
interaction with the help of a geometrical capacitance C 
the incremental charge is simply related to the poten- 
tial variations in the conductor 5U and the gate 5V g via 
5Q = C(SU — 5V g ). On the other hand, the charge SQ is 
also the sum of the injected charges due to the variation 
of the contact potentials (keeping the internal electro- 
static potential fixed) and the induced charges which are 
generated by the Coulomb interaction. Thus 



SQ = C(SU - SV g ) = i [q p N p (p)Sfi p + q h N h (p)S f i p 

+ q h N h {h)5 l l h + q p N p {h)8 i l h 

- q p N p {p)q p 8U p - q p N p {h)q h 5U h 

- q h N h (p)q p 6U p - q h N h (h)q h 5U h ] (5) 

The first four terms in Eq. (||) determine the charge 
injected into the conductor as a consequence of the varia- 
tion of the contact potentials. Similarly to the conceptual 
separation of the electrostatic potential, we have here as- 
sumed that each contact has separate electrochemical po- 
tentials fj, p and fjih for electrons and holes Q. In reality 
Hp = —fih = A* = eV'. The second four terms in Eq. (|5|) 
determine the induced charges. In the absence of super- 
conductivity Eq. (0) reduces to the "Poisson equation" 
for charge variations in a normal conductor 00]: We 
have N p (h) = N h {p) = and N = N p {p) = N h {h). Go- 
ing back to Eq. (|B|) we can solve it for SU and can use this 
solution to determine the electrochemical capacitance of 
the hybrid structure vis-a-vis the gate. We find 



r<NS _ &Q 
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5Q Ce 2 N^ 
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with a total density of states 

Nz = \[N v {p) - N p (h) + N h (h) - N h (p)}. 



(6) 



(J) 



It is instructive to examine C,, for an ideal ballistic 
wire. For a perfect N — S structure, C p is zero, since ev- 
ery charge incident upon the conductor is perfectly An- 
dreev reflected and hence the net accumulated charge is 
zero. If the superconductor is driven into the normal 
state, the capacitance is C% = Ce 2 N/(C + e 2 N) 

which for a ballistic wire reduces to — C since typi- 
cally e 2 /C » 1/N. Thus C^ s for the hybrid structure 
can differ dramatically from . 

If we now wish to consider fluctuations in the charge 
then we have to consider the Poisson equation for the 
fluctuating charges H,||. For the case that the voltages 
are held fixed (zero-impedance external circuit) this leads 
to an operator equation, 



Q = CU = eJV - e 2 N^U 



(8) 



where Af is the operator of the bare charge fluctuations 
given by Eq. ^ and the last term in Eq. (||) describes 
the screening of the bare charge fluctuations. Solving 
Eq. (||) for U we can express the charge fluctuations 
Q in terms of the bare charge fluctuations {I. For the 
fluctuation spectrum of the charge this gives 



S QQ M = (l/2)^iV E 2 ^]r J dEF 7ah (E,hu) 



7<5 a(3 



; Tr[AT laS , (E, E + M^U, (E, E + Tiuj)] 



(9) 
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where 



N laS ,{E,E') = J2sgn(q^fd 3 rUl 5 (^,r;E 7 E'), (10) 



F^ Sp (E > E') = f la (E)[l-fs,(E')} 
+fs (E')[l-f la (E)}. 



(11) 



In the above the sum is over all normal contacts 7<5 and 
degrees of freedom a/3. The Fermi functions f la (E) 
are defined such that / 7 = fo(E p — u 7 ) and f lh = 
fo(Eh + M7) where E a is the energy of a particle (hole) 
(a = p(h)) in reservoir 7, which is at a chemical poten- 
tial fij, fo(E) is the Fermi function at the condensate 
chemical potential of the superconducting leads (/io)- 

Now we evaluate Eq. (||) at equilibrium and zero tem- 
perature, to leading order in hu> <C |A|. We find 



S QQ (u) = 2ClR q h\w\ 
with a charge relaxation resistance 



R„ 



4e 2 [7V S ] 2 



(12) 



(13) 



In the zero frequency limit, and at a temperature kT -C 
|A| the charge fluctuations are given by Sqq(uj) — 
2C 2 R q kT with R q as given by Eq. (|l|). If we bring 
the hybrid structure into a non-equilibrium state by ap- 
plying a bias e\V\ ^> kT between the normal reservoir 
and the superconductor, we find at zero temperature in 
the low frequency limit to leading order in e\V\ for the 
system shown in Fig. [l], 

S qq (cj) = 2€%Rve\V\ (14) 

with a non-equilibrium resistance 

h TiM,H+^A) 



4e 2 



[N T 



(15) 



The resistance Rv is a consequence of the charge fluctu- 
ations which arise due to the shot noise of the transport 
state. The fluctuation spectra of the current at the gate 
G in Fig. (Q) are related to the charge fluctuation spectra 
Eqs. (|f) and (H3) via S n (u) = lu 2 S qq {u)_ 

To proceed further we follow Beenakker || and express 
the electron-hole scattering matrix of the NS-structure 
in terms of the scattering matrix elements of the normal 
structure. For simplicity we restrict our considerations 
here to normal structures which are symmetric with re- 
spect to left and right going carriers. The scattering ma- 
trix of the normal conductor is then given by a reflection 
matrix r and a transmission matrix t only. In terms of r 
and t the electron-hole scattering matrix elements are 

s lplp (E) = r p (E) + a 2 t p {E)r h {-E)M p {E)t p {E) (16) 

s lhlp (E) = ae-^t h (-E)M p (E)t p (E) (17) 

s lhlh {E) = r h {E)+a 2 t h (E)r p {~E)M h (E)t h {-E) (18) 

s lplh (E) = ae^t h (E)M h (E)t p (-E) (19) 



where a = exp[— iarccos(E/A)] and 

M p (E) = [l - a 2 r p (E)r h (-E)] 
M h {E) = [l-a 2 r h (E)r p (-E)] 
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(20) 
(21) 



In Eqs. (|16|-21 ) we have given an index p, h to r and t 
to indicate whether a carrier or a hole is scattered since 
we have to distinguish derivatives with respect to the 
electron and hole potentials. Throughout the rest of the 
work presented here we will be making use of particle- 
hole symmetry (ie. E — > s p = s* h ). 

Now we wish to calculate the properties of the sys- 
tem shown in Fig. 1. We choose to model a quantum 
point contact (QPC), using a saddle point potential JO] 
V(x, y) = Vo + ^muj 2 y 2 — ^muj 2 x 2 and as in Ref. || we 
have evaluated the density of states semi-classically. We 
also take the WKB limit allowing us to make the fol- 



lowing transformation 



d/dE x . In 



i!)" 1 dq x U x (v) 

the above transformation we have been careful to keep 
the index A included, such that we can distinguish be- 
tween particles and holes. We also note that we can write 
ds/dE x = {ds/d9 x ){dO x /dE x ) where dO x /dE x is the par- 
ticle (hole) (A = p(h)) density of states divided by n in 
the conductor in region Q. Making use of particle- hole 
symmetry we have dO/dE = d6 p /dEP = d9 h /dE h . 
Using Eqs. ( |Il||l6f]l^ ) we now calculate R q and find 
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NS 



h E, 



16-R„+-R„(4fl„ + l) 2 r^ 
(l + fl„) 4 



\ dE J 



2e 2 



(22) 



as compared P|l3| to 



R. 



N _ h En ( dE ) 

^ [£ n (^)f 



(23) 



for a QPC attached to two normal leads, where T n and 
R n = 1 — T n are the transmission and reflection proba- 
bilities of the nth quantum channel and 9 n is the phase 
accumalated by a carrier in the nth channel during trans- 
mission through the QPC. In Fig. | we plot R^ s and 
R^ (the charge relaxation resistance for a QPC attached 
to two normal leads). For the parameters we have cho- 
sen, we see that the difference between and R^ s is 
roughly a factor 2. This arises from the difference be- 
tween the contact resistance of a mesoscopic conductor 
in an N-S system compared to an N-N system |l4| ]. 

Having examined R^ s and R q we now proceed to ex- 
amine Ry , finding 



r>NS _ 
K,/ — 



h E, 



(l + fi„) 4 V dE . 



2e 2 



E f 



d6 n \ 
dE I 



(24) 



3 



0.02 



0.01 



0.00 





f 

/ 


JL , ' 


1 


J 



0.0 



1.0 



E _V o 2.0 



3.0 



6.0 

4.0 
G 
2.0 

0.0 



1e-03 



FIG. 2. R^ s (thick solid line) and R% (thick dashed line), 
in units of h/e 2 . G NS (thin solid line) and G N (thin dashed 
line) is the conductance, in units of e 2 /h, as a function of 
Ef — Vo in units of hu> y , for a saddle QPC with lu v /lu x = 2. 



as compared M to 



T>N 



h E 



n 4R„T r , 
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(25) 

e En ( He)] 
for a QPC attached to two normal leads. In Fig. || we 
see both Ry and Ry as a function of Ep — Vq. At the 
channel threshold R n — T n — 1/2 the non-equilibrium re- 
sistance Ry of the hybrid structure tends to (9/4) (h/e 2 ) 
whereas the normal state Ry s is completely suppressed 
at the channel threshold due to the divergence of the 
semiclassical density of states. 

In this work we have taken the first steps in considering 
frequency-dependent charge fluctuations and noise cal- 
culations in mesoscopic normal superconducting hybrid 
systems. To do this we have considered structures where 
the condensate chemical potentials of any superconduct- 
ing lead remains the same and gates only see the normal 
regions of the conductor. This allows us to consider os- 
cillating voltages in the normal reservoirs and calculate 
the systems response. We have focused on the i?C-time 
of such structures. The i?C-time is a fundamental dy- 
namical quantity of electrical conductors. Elsewhere we 
have shown that for normal mesoscopic conductors the 
i?C-time also determines the dephasing rate of Coulomb 
coupled structures jL{|. We expect that such a relation 
also holds for hybrid systems. As an example we con- 
sidered the most simple non-trivial example, a QPC at- 
tached to one normal lead and one superconducting lead. 
We have given a firm prediction that the i?C-time and 
the fluctuations induced into a gate differ markedly for an 
NS structure and a normal conductor. As usual, we as- 
sume perfect Andreev reflection at the superconducting 
interface, whereas at real interfaces only partial Andreev 
reflection might occur. Currents induced into gates can 
be measured |l6| . We believe that it is similarly possible 
to measure the current fluctuations at a gate. 
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FIG. 3. Ry s (solid line) and R$ (dashed line), in units 
of h/e 2 , as a function of Ef — Vo in units of hui y , for a saddle 
QPC with cjy/ujx = 2. 
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